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ABSTRACT 

The full stellar population of NGC 6334, one of the most spectacular regions of 
massive star formation in the nearby Galaxy, have not been well-sampled in past studies. 
We analyze here a mosaic of two Chandra X-ray Observatory images of the region using 
sensitive data analysis methods, giving a list of 1607 faint X-ray sources with arcsecond 
positions and approximate line-of-sight absorption. About 95% of these are expected to 
be cluster members, most lower mass pre-main sequence stars. Extrapolating to low X- 
ray levels, the total stellar population is estimated to be 20 — 30, 000 pre-main sequence 
stars. The X-ray sources show a complicated spatial pattern with ~10 distinct star 
clusters. The heavily-obscured clusters are mostly associated with previously known 
far-infrared sources and radio HII regions. The lightly-obscured clusters are mostly 
newly identified in the X-ray images. Dozens of likely OB stars are found, both in 
clusters and dispersed throughout the region, suggesting that star formation in the 
complex has proceeded over millions of years. A number of extraordinarily heavily 
absorbed X-ray sources are associated with the active regions of star formation. 

Subject headings: ISM: clouds - ISM: individual (NGC 6334) - stars: early type - stars: 
formation - stars: pre-main sequence - X-Rays: stars 



1. INTRODUCTION 

The NGC 6634 (= Sharpless 8 = RCW 127) complex of molecular clouds and HII regions is an 



unus ually rich region of massive st ar formation in the nearby Galaxy (see review by lPersi Tapia 



2008). At a distance of 1.7 kpc (INeckell Il978l ) and Galactic coordinates (351.2+0.7), it exhibits 
many phenomena associated with embe dded OB stars: luminous far-in frared sources denotin g 
embedded clusters (jMcBreen et al.l 1 19831): radio continuum H I I regions (jRodriguez et al.l [l982j) ; 
and complexes of molecular masers (jMoran fc Rodriguez! Il980i ) . Early studies argued that the 
star formation started in the center of the region, where the stellar cluster appears dynamically 
relaxed, and proceeded outward where the embedded protostars are seen. Recent high-resolution 
millimeter-band maps of the dust distribution reveal hundreds of molecular clumps distributed 
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along a narrow ridge along the Galactic Plane (jMunoz et al.l 120071 ; iMatthews et al.l |2008| ) . The 
total molecular mass is estimated to be ~ 2 x 10 4 Mq with typical column density ~ 5 x 10 22 cm -2 
along the ridge. The total bolometric luminosity is ~ 2 x 10 6 L ( Loughran et al. 198d ) indicating 
that rich embedded stellar clusters with massive OB stars are present. Optical and infrared maps 
are dominated by bright and spatially complex nebular emission from HII regions and heated dust. 

While deeply embedded massive stars in the region have been well-studied, less is known 
about the low mass pre-main sequence population associated with the massive stars. JHK source 
lists are incomplete due to heavy obscuration and H II region nebulosity. A wide field survey 
of 11 S, K < 13.5 sourc es finds hundreds of embedded low mass stars around the high mass 
regions (IStraw etaL[|l989|) and port ions of the cloud have been surveyed to fainter levels in JHK 



(jTapia et al 



1996 



Persi et al.ll2000l ). At both bright and faint magnitudes, near- infrared samples 



are severely contaminated by unrelated Galactic field stars. 

High-resolution X-ray images of star forming molecular clouds from the Chandra X-ray Obser- 
vatory can give rich samples of p re-main sequence sta rs which are complementary to those found 
in infrared surveys (see review bv lFeigelson et al.ll2007l ). The X-ray band reveals all types of young 
stars. OB stars emit X-rays from their stellar winds on several scales: small instabilities in the 
wind acceleration region close to the stellar surface; magnetically channeled wind shocks; colliding 
wind shocks of close binaries; and shocked winds on parsec scales suffusing large H II regions (e.g . 



Gohen et alihoosl : iTownslev etldlhoO.i Isdiulz et al 



2003 



Gagne et al 



2005 



Skinner et al 



2005 



Sana et al.1 120061 ) . Surveys of the Orion Nebula Cluster and Taurus clouds show that all lower 
mass stars emit X-rays during magnetic reconnection even ts near the stellar surface similar to, bu t 



far more powerful and frequent than, flares on the Sun (jGetman et al 



2005 



Giidel et al.ll2007f ). 



Studies with the high-resolution Chandra X-ray Observatory are particularly effective in produc- 
ing stellar samples with several benefits: low contamination by unrelated older stars; penetration 
deep into the obscuring cloud; sampling of rich populations of lower mass stars; (sub)arcsecond 
positions; and selection that is mostly unbiased with respect to the presence or absence of proto- 
planetary disks. Chandra studies of rich stellar c lusters in obscur ed star form ing regions aroun d 
1.5 - 2 kp c distance include the Trifid Nebula feho et al.1 boodl RCW 38 dwolk et alJhoOfih . 
Cvg OB2 (Ulbacete Colombo et al.l [2007l). M 17 (|Broos et al.l boojl). NGC 6357 ne ar NGC 6334 
(|Wang et al.ll2007l ). and the Bosette Nebula and Molecular Cloud (j\Vang et al.l [20081 . bood ). 



Two Chandra images of NGC 6334 are presented bv lEzoe et al.l (|2006l ). Their study emphasizes 
the hard X-ray emission from individual stars and possibly from diffuse plasma suffusing the HII 
regions. They extract a sample of 792 X-ray sources, but do not provide a source list to facilitate 
further study . Prior to the Chandr a observation, low resolution X-ray images gave only a glimpse 
of structure (ISekimoto et al.1 120001) . Ultra-h ard X-ray emission attributed to background active 
galaxy has also been seen (jBykov et al.l l2006). 



We report here a reanalysis of the two Chandra exposures of NGC 6334 examined bv lEzoe et al 



(J2006J) using advanced data analysis techniques. Section [2] describes the data analysis methodology, 
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Ogives the astronomical results, and §H presents conclusions. We find 1607 faint X-ray sources of 
which only ~5% are likely problem sources or non-cluster contaminants. We give accurate positions 
and estimated absorption measurements for each star. From this dataset, we obtain an improved 
view of the clustering in the region: 10 X-ray clusters are identified, including several previously 
unreported lightly obscured star clusters. A considerable number of massive OB stars, some newly 
found, and deeply embedded protostars are detected. These may include extremely young Class 
massive protostars. 



X-RAY DATA ANALYSIS 



2.1. Chandra Observations 



The X-ray observations discussed here were previously presented bv lEzoe et al.l (|2006l ). NGC 
6334 was observed in two overlapping exposures of the Imaging Array of the Advanced CCD 
Imaging Spectrometer (ACIS-I), each 40 ks in duration and covering a ~ 17' x 17' field of view. 



The Chandr a X-ray Observator y satellite and ACIS instrument are described by IWeisskopf et al 



(|2002l ) and lGarmire et al.l (|2003l ). The Chandra mirrors give an excellent subarcsecond point spread 
function (PSF) on-axis, but the PSF degrades substantially with off-axis distance. The field centers 
are (a, 5) = (17 h 20 m 01.0 s , -35°56'07") and (17 ft 20 m 54.0 s , -35°47'04") with a satellite roll angle of 
269°. The exposure times are 40.2 ks and 40.1 ks, respectively. The observation was made in the 
standard Timed Exposure, Very Faint mode, with 3.2 s integration time and 5 pixel x 5 pixel event 
islands. Pre- Chandra X-ray obse rvations with the low-resolution ASCA satellite are presented by 
Sekimoto et al.l (|2000l ). They and lEzoe et al.l (|2006l ) examine the spatially- integrated X-ray spectra 
of different components of NGC 6334 which we do not address here. 



2.2. X-ray Source Identification 



Our data reduction procedure follows those of iTownsley et al.1 (120061 ) and iBroos et al.1 (120071 ) 
with a few refinements. Level 1 event l ists a re first processed to recover an improved Level 2 
event list as described by lTownsley et al.l (120031 . Appendix B). This includes filtering of bad events, 
removal of the ±0.25" positional randomization, and the correction of charge transfer inefficiency. 
Cosmic-ray afterglows are removed by searching for multiple events arriving on the same detector 
pixel within a few CCD readout frames with decreasing energ ies. Sub-pixel po sitioning of events 



with non-zero CCD 'grades' are refined using the algorithm of (jMori et al.ll200ll ). Matching bright 



X-ray sources in the central region of the detectors to stars in the 2MASS Point Source Catalog, 
a 0.35" correction was applied to astrometrically align the fields to the Hipparcos reference frame. 
Figure [H shows the resulting merged ACIS fields at reduced resolution, and Figure [2] shows two 
expanded views. 
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It is challenging to reliably identify faint sources in overlapping, crowded fields with a spatially 
varying PSF, and to establish their statistical significance. We adopt a procedure where a super- 
set of possible sources is first constructed by a variety of means, following by iterative pruning of 
low-signi ficance sources. An i nitial list was obtained by running the wavdetect wavelet detection al- 
gorithm ( Freeman et al.ll2002l ) with a low threshold of 10~ 5 . Other potential sources are added from 
visual inspection of the image s and from a peak-finding algorithm applied to a maximum-likelihoo d 
reconstruction of the images (|Lucylll974h using local PSFs, as described by lTownsley et al.l (|2006l ). 
The Lucy-Richardson reconstruction is applied to the crowded central arcminute of the NGC 6334 
fields with 0.1" pixels. These steps are particularly important to identify closely-spaced sources 
which are merged by the wavelet algorithm. These procedures produced a superset of 2511 candi- 
date sources. A few sources missed by the algorithm found by visual examination were added. 

Photons are then extracted in polygonal regions closely matched to the local PSFs for each 
candidate source using the IDL-based program ACIS Extract. The number of extracted events in 
the total Chandra 0.5 — 8 keV band is compared to the local source-free background level. The 
statistic Pb measures the probability that the extracted eve nts would be found from random Poisson 
fluctuations in the local background (see the Appendix of IWeisskopf et al.l 120071 ) . This use of Pb 
as a source existence criterion has the advantage of being independent of the poorly understood 
sensitivity limits of the wavdetect and maximum-likelihood restoration processes in the presence of 
crowding and spatially varying PSFs. 

For our analysis here of the NGC 6634 fields, we choose a significance threshold of source 
existence Pb < 0.001 (logPe < —3.0) for source existence. A list of candidate sources with 
Pb < 0.01 was constructed to obtain measurements of local source-free background levels for the 
Poisson probability calculation. For each source, Pb values are calculated in three X-ray bands 
(0.5 — 2 keV, 2 — 8 keV and 0.5 — 8 keV) and, in the region overlapping the two ACIS fields, in both 
fields, and the lowest value is kept. After four iterations of this procedure in ACIS Extract, where 
candidate source significances are recalculated based on revised background levels, we emerge with 
1607 sources with Pb < 0.001. These sources are listed in Table CD Note that we obtain twice 
as many sources as the 792 sources obtained by Ezoe et al. (j2006l ) using traditional data analysis 
methods. 



Table [T] has a format similar to the source tables in iGetman et al.l (J2005J) and iTownslev et al 



(|2006l ). Column (1) gives the ACIS running sequence number, and column (2) provides the IAU 
designation. Columns (3) — (6) show R.A. and decl. in degrees, positional error in arcseconds, and 
off-axis angle in arcminutes. Columns (7) — (11) give net counts in the full band, the associated 1 a 



1 ACIS Extract (AE) and its associated Tools for ACIS Review and Analysis (TARA), developed by P. Broos and 
colleagues by the ACIS Team at Penn State, are publicly available at http://www.astro.psu.edu/xray/docs/TARA 
AE is a versatile script written in IDL that performs source extraction, facilitates fitting X-ray spectra, creates light 
curves, and computes a wide variety of statistical properties of the sources. This analysis was performed using AE 
version 3.79. 
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equivalent uncertainty based on Poisson statistics, background counts in the extraction region, net 
counts in the hard band, and the fraction of the PSF used for source extraction. PSF fraction values 
below ~ 0.9 indicate source crowding. Column (12) presents the Pb probability of photometric 
significance when they lie in the —5 < log Pb < —3 range. 

Column (13) of Table [1] gives the sources median energy E me d in the full spectral band after 
background subtraction. The median energy is a reliable surrogate for interstellar gas absorption. 
For E me d > 1.7 keV, an empirical relationship with the gaseous column density lo giVff = 21.22 + 



0A4E me d cm 2 is found for pre-main sequence stars in the Orion Nebula Cluster (jFeigelson et al 



20051 . see their Figure 8). This relationship is also found in a large population of M 17 stars 
(Getman et al., in preparation) and we assume it applies for NGC 6334 stars here. For normal 
interstellar gas-to-dust ratios, the column densit y can be convert ed to a visual absorption with 
the traditional relation Ay = Nn/2 x 10 21 mag ( Gorensteinlll975l ) o r the revised relati on derived 
from X-ray absorption measurements, Ay = Nh/1.Q x 10 21 mag ( Vuong et al.1 120031 ). We use 
these calibrations between the instrumentally dependent E me( i and physical absorption measures 
repeatedly in the science sections below. 



2.3. Stellar Counterparts 

Unfortunately, no wide-field, sensitive, and high-resolution near-infrared image has been pub- 
lished covering the ACIS fields of view in NGC 6334. We thus use the 2MASS Point Source Catalog 
to find stellar counterparts to the X-ray sources. Column (14) in Table Q] gives the K s magnitude 
of 2MASS stars lying within I" of the X-ray source positions. About half of the Chandra sources 
have a 2MASS star within 1", but the majority of these have inaccurate photometry due to nebular 
contamination, star crowding, or sensitivity limits. 

Figure [3] shows JHK color-color and color-magnitude diagrams for the 2MASS stars with 
photometric errors better than ±0.1 mag in all three bands. The results are similar to diagrams 
obtained from other C/iandra-selected populations in clusters with distanc e s aro und 1.5 — 2 kpc 



such as M17 and the Rosette Nebula (IBroos et al.1 120071 : IWang et al.ll2008l . [2009J). The majority 



of X-ray selected stars do not exhibit i^ s -band excesses attributed to the heated inner edges of 
protoplanetary disks; most would be classified as Class III pre-main sequence stars based on near- 
infrared photometry. The 2MASS photometric diagrams in Figure [3] show a wide range < Ay < 
25 mag. Noting that the X-ray hardness measure E me d exceeds 3.0 keV (log Njj > 22.5 cm~ 2 or 
Ay > 20 mag) for many of the X-ray sources without 2MASS counterparts, stellar counterparts 
for many X-ray sources will be too heavily absorbed to be present in the 2MASS catalog. Our 
experience with M17 and Rosette where deeper and higher-resolution near-infrared imaging was 
available also indicates that the 2MASS catalog includes only a fraction of the true stellar population 
represented by the Chandra sources biased towards those with higher bolometric luminosities, lower 
obscuration, and locations away from bright nebular emission. Due to these limitations of the 
2MASS counterparts, we do not pursue scientific interpretation of these near-infrared sources (e.g., 
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disk fractions, Initial Mass Function). 

The final column of Table[T]gives possible associations between X-ray sources and objects from 
published objects in NGC 6334. These were obtained from the SIMBAD bibliographic database 
using a ±3" pos itional agreement . The most common possible counterparts are from the shallow 
JHK survey of IStraw fc Hylandl (119891 ) with labels FIR-I to FIR-V corresponding to the major 
far-infrared clusters. We emphasize that some of these associations may not be correct due to 
positional inaccuracies or offsets between the X-ray and long-wavelength bands. 



2.4. Completeness and Contaminants 



Our source detection method is complete down to a statistical probability limit which can be 
converted to an ACIS count limit and, with considerable uncertainty, an X-ray flux limit. The 
faintest on-axis sources in Table Q] have 3 — 4 counts. Based on comparison with ot her pre-main se- 



quen ce studies of pre-main sequence populations (e.g., the sensitivity equation (7) in lFeigelson et al 
2002 ). this corresponds to an absorption-corrected luminosity in the 0.5 — 8 keV of Lt c ~ 4 x 10 29 



erg s _1 for absorption around log Nh ~ 21 — 22 cm -2 . This limit is quite approximate, as it depends 
on the intrinsic spectrum and cloud absorption; the limit will be a factor of ~ 2 lower for lightly- 
absorbed sources, several times higher for sources with log Nh ~ 23 cm -2 , and roughly ~ 100 times 
higher for sources with log Njj ~ 24 cnr 2 . Significant instrumental spatial variations in sensitivity 
are also present due to off-axis telescope vignetting and degraded point spread function in each 
field and to the overlap region between the two exposures. While it is thus difficult to establish a 
single X-ray completeness limit for the NGC 6334 observations for these reasons, we are confident 
the most of the observed region is surveyed to a sensitivity better than 1 x 10 30 erg s _1 except for 
sources suffering extremely high absorption. 



Based on the well-established correlation between X-ray luminosity and mass (jTelleschi et al 



20071 ). this sensitivity limit is sufficient to c apture most of the cluster members with masses above 
1 M Q (except for a minority of A and B stars. IStelzer et al. 2006), and some stars between 0.3—1 M . 
A survey complete down to ~ 1 M Q detects about 5% of t he stellar popu lation (including young 
brown dwarfs) assuming a standard Initial Mass Function (IChabrierll2003l ). We can thus multiply 
the number of X-ray source by a factor of ~ 10 — 20 to estimate the total population of each cluster. 
This is only an approximation, as the X-ray sensitivity varies with absorption across the field and 
in the radial direction. 

Our method of locating Chandra sources is considerably more sensitive than that used by 
Ezoe et al.1 (|2006l ). but it entails several sources of error. Careful examination of the image suggests 
that several sources located in the wings of strong sources may be artifacts of the processing. 
We cannot immediately exclude these as spurious because studies of other young stellar clusters 
establish that close mult iple X-ray source s can appear around massive stars (see HD 46150 in 
the Rosette's NGC 2244, IWang et al.ll2008l ). Some of the faintest sources are also likely not real, 
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particularly in the densest cluster regions. Recalling that ~ 10 — 20 very X-ray faint stars should 
be present for each detected source, many low mass clusters members are expected to produce 
1 — 3 counts which are viewed as background by ACIS Extract. I n addition, t r uly d iffuse plasma 
emission from shocked O star winds may present, as discussed by lEzoe et al.l (|2006l ) . The result 
of the elevated background is to reduce the statistical significance of real sources and to cause the 
values of log Pb to fluctuate with each iteration of A CIS Extract as different candidate sources are 
rejected and their photons enter the background. This effect generally leads to a more conservative 
list of faint sources than an environment without emission from diffuse plasma or unresolved stellar 
populations. Nonetheless, very weak sources in the most crowded regions have an increased chance 
to be spurious. Altogether we believe that no more than a few percent of the sources in Table Q] are 
spurious. Improved methods for pruning the source list in crowded regions are under development 
(Broos et al., in preparation). 

A small fraction of the reliable X-ray sources are not members of the NGC 6334 population. 
Detailed modeling of contaminant populations in the foreground or ba ckground of Ga l actic Plane 
molecular clouds at distances around 1 — 2 kpc have been made by iGetman et al.l (120061 ) and 



Wang et al.l (|20071 ). Scaling these results to the exposures, survey area and distance of NGC 6334, 
we estimate that 25 — 50 sources are extragalactic sources (mostly quasars), 15 — 30 are foreground 
field stars, and < 10 are background field stars. Thus, 2 — 5% of the 1607 sources in Table [T] are 
probably contaminants. However, these sources must be randomly distributed across the field and 
their contribution to the stellar clusters discussed here will be negligible. 



2.5. Defining X-ray clusters 



Cursory examination of the ACIS images in Figures [T][2] shows that the spatial distribution 
of X-ray sources is very inhomogeneous with distinct strong co ncentrations that represent young 
stellar clusters. This was previous l y foun d by lEzoe et al.1 (120061 ): see their Figures 2-3. Adopting 
the procedure used by IWang et al.l (|2009l ) , we locate clusters on a map of the X-ray source density 
obtained by smoothing the locations of the 1607 sources with a Gaussian kernel of width a = 45" 
(0.4 pc). Here we ignore the few sources that are likely non-members of the cluster, and consider the 
X-ray source surface density to represent the stellar surface density of the NGC 6334 population. 
Note that these are not maps of smoothed X-ray emission; rather, the X-r ay stars are viewed 
individual stars tracing the spatial structure of the full young stellar population. IWang et al.1 ([2009 ) 
show that the existence of structures smaller than several arcminutes is not strongly affected by 
inhomogeneities in the X-ray sensitivity across the mosaic due to reduced sensitivity far off-axis of 
each pointing and increase sensitivity in the overlap region between adjacent images. Nonetheless, 
these variations in sensitivity reduce the accuracy of our estimates of cluster richness. 

The resulting smoothed stellar surface density map is shown in FigureHJfor the total 0.5 — 8 keV 
Chandra band, and in Figure [5] for a 'soft' and and 'hard'. Here we choose to divide the sample 
into soft sources with E me d < 2.5 keV and hard sources with E me d > 2.5 keV; the boundary 
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approximately corresponds to column density log Nh ~ 22.3 cm -2 and visual absorption Ay ~ 10 
mag. Hereafter, we refer to these subsamples as the 'lightly obscured' and 'heavily obscured' X-ray 
stars and clusters. 

We have chosen to define clusters as localized peaks in X-ray source surface density which are 
at least 10-fold above the background level. This corresponds to the yellow contour in Figure [H 
where dark blue represents the background level. Eight or nine clusters can be discerned in this 
total band map. However, examination of the smoothed soft and hard band stellar surface density 
maps shows dramatic differences; most clusters are prominent only in one band (Figure [5]). From 
examination of the three smoothed maps, we emerge with a list of ten clusters listed in Table [2] 
and labeled in Figure [5j The table columns give our designations XA through XJ, central location 
(accurate to ±20"), indicator of X-ray hardness, approximate extent based on the yellow contour in 
Figure HI and the number of X-ray sources found inside this contour for each X-ray cluster. Recall 
that 1' corresponds to about 0.5 pc at a distance of 1.7 kpc. 

The procedure used here to define X-ray selected clusters based on peaks in a smoothed dis- 
tribution is not unique. Other statistical techniques can be used: sm oothing using different kernel 
widths; nearest-neighborhood methods (jRoman-Zuniga et al.ll2008l ). single-linkage agglomerative 
clustering (also called the 'friends-of-friends' algorithm), maximum-likelihood mixture models as- 
suming Gaussian cluster shapes, and so forth. E ach of these metho ds has advantages and disad- 
vantages, and none gives a unique 'correct' result. IWang et al.1 (|2009l ) compares our procedure to a 
10-th nearest neighbor algorithm in a Chandra mosaic of the Rosette Molecular Cloud, and shows 
that the results are very similar. For any clustering method, the results may not be physically 
meaningful. For example, it is possible that the four clusters XA-XD represent a single physical 
structure that appears divided by inh omogeneous ob scuration. The confusing effects of variable 
extinction in NGC 6334 was raised by iLindsavl (|1955l ) and others. Our trials in constructing the 
cluster list with different smoothing and associations between hard and soft structures gave 8 — 12 
clusters. We conclude that the division of the X-ray sources into ten distinct clusters is a reasonable, 
but far from unique, result. 



2.6. Comparison to Long- Wavelength Studies of NGC 6334 



Figure [6] superposes the smoothed X-ray source distributions on wide-field maps at radio 
and far-infrared w avelengths. The top panels show contours from the Very Large Array 18 cm 
co ntinuum map of Sarma et al.l (|2000l ) with labels for the NGC 6334 A-F HII regions identified 



by iRodriguez et al 



Carral et al 



(|1982l ). Higher resolution radio maps of some co mponents are p r esente d by 
(2002). The bottom panels show a map at 71 /iin by iLoughran et al.l (|1986l ) ob- 
tained with a balloon-borne telescope with their designations of far-infrared sources I-V. Somewhat 
higher resoluti o n far -infrared maps from HIRES analysis of the IRAS satellite data appear in 



Kraemer et al 
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Examination of these map overlays shows that there is a poor correspondence between lightly 
obscured X-ray clusters and radio/infrared structures, but a good correspondence when heavily 
obscured X-ray clusters are considered. Thus, XE appears associated with far-infrared IV and 
radio A, XH is associated with II and D, and XI is associated with far-infrared I and radio E 
and F. The lightly obscured clusters XA-XD, XF, XG, and XJ have little association with long 
wavelength structures, although some of the stars we include with XC-XD appear embedded and 
may be associated with far-infrared V and the southwestern radio shell. 

Figure [7] shows the X -ray clusters superpos ed on the mid-infrared emission observed by the 



Spitzer GLIMPSE survey ([Benjamin et al.ll2003l ). Some association can be seen with the heavily 



obscured clusters: XE lies in a confused region of bright and dark infrared nebulosities; XH is 
straddled by large bubbles emerging from the ridge line; and XI coincides with an isolated bright 
region within a larger dark cloud. There is also a hint that the massive stars in lightly obscured 
X-ray clusters (particularly XA, XC, XD and XG) have produced bowl-like cavities in the cloud 
material traced by the Spitzer map. But the relationship between the X-ray maps (which mostly 
represent lower mass stars) and the mid-infrared maps (which mostly represent heated dust and 
PAH molecules) is indirect. We do not significantly use the GLIMPSE map in the science analysis 
of this study, as we focus on the stars rather than the interstellar material. 

We proceed with a more detailed examination of individual stars and star groups in various 
components of NGC 6334, proceeding west to east. 



STARS AND STELLAR CLUSTERS IN THE CHANDRA POPULATION 
3.1. Embedded Stars in NGC 6334 V 



A group of faint embedded X-ray sources around (a, 5)=(17 /l 19 m 58 s , -35°58') lying behind 
the lightly absorbed cluster XD can be seen in the hard band image. Its surface density was not 
high enough to be classified as an X-ray cluster. The group includes sources #218, 219, 226, 228, 
234, 236, 248, 249, 252, 266, 270, 282 and 294. Several of the sources exhibit extraordinarily 
high X-ray absorptions with 4 < E me d < 7 keV, corresponding to 23.0 < log Njj < 24.3 cm -2 
and 50 < Ay < 1000 mag based on the calibration discussed in §2.21 In nearby star formation 
regions, only a handful of young stellar objects with dense localized envelopes ha ve comparable X 



ray absorptions including the Class 0/1 sources CrA IRS 7 and IC1396N BIMA 2 (jHamaguchi et al 



2005 



Getman et al.ll2007l ). 



These X-ray sources lie on the southeast side of the NGC 6334 V far-infrared region and 
surround the embedd ed B-type protostars p roducing masers and the infrared emission nebulae 
IRN V-l and IRN V-2 lHashimoto et al.l (J2008J). These massive protostars are too heavily obscured 
to be seen in the near-infrared, consistent with the extreme absorption detected in the nearby X-ray 
sources which probably represent the most X-ray luminous low mass members of the associated 
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cluster. This grouping bears some resemblance to the cluster of ~ 20 X-ray emitting low mass stars 
associated with the cluster of B-t ype protostars in the Kleinman-Low Nebula behind the Orion 
Nebula with logN H ~ 23.0 cm" 2 (|Grosso et alJbood ). 



A sparse but compact (diameter < 0.1 pc) group of faint, moderately absorbed X-ray sources 
appears to be associated with the bright Ha emission line star SS 303 (K=8.57) and the infrared 
source NGC 6334 V IRS 1. These include # 302, 305, 306, 310, 311 (= SS 303) and 320 (= 
NGC 6334 V IRS1). iJackson Kraemerl (|1999l ) suggest that SS 303 ionizes the the southeastern 
radio shell, G351.20+0.70. Another possible ionizing source is the bright star #210 = FIR-V 51 
(K=7.76), although it lies off-center within the radio shell. 



3.2. Lightly Obscured Clusters XA-XD 

The XA, XB, XC and XD clusters together constitute the largest and richest complex seen 
in the NGC 6334 region in the X-ray band. It is unclear whether they are structurally distinct 
structures or part of a single larger cluster subject to inhomogeneous absorption. The latter in- 
terpretation is supported by the smoother appearance of the source distribution in the hard band 
than in the soft band (Figure [5]). These X-ray sources are lightly absorbed with E < 2.5 keV, and 
relatively few are present in the hard X-ray maps. This obscuration corresponds to log Nh < 22.3 
cm -2 and Ay % 10 mag. 

These X-ray structures have few published counterparts in other wavelengths, perhaps because 
they are older clusters and no longer embedded in dense cloud cores. The mid-infrared source 
IRAS 17175-3554 corresponds closely with XC, although we suspect is lies behind the cluster in 
the cloud interior. C luster XD can be seen in the near-infrared H2 and L band images obtained by 



Burton et al.l (|2000l . see their Figure 5) where it was associated (perhaps incorrectly) with the em- 



bedded far-infrared source NGC 6334-V. Many of the X-ray sources with 2MASS c ounterparts with 



11 < K < 14 are likely low mass cluster members, and some are FIR-V stars from lStraw &: Hvland 



(|1989l ). A number of stars bright in both the X-ray and K s bands are likely high mass cluster 
members. These include # 210 (= FIR-V 51, see i £U|) . 354 (=FIR-V 10), 405 (= FIR-V 102), 415 
(= FIR-V 107), 505, 564, and 586. Source #405 = FIR-V 102 (K s = 8.81) lies in a compact group 
(diameter ~ 0.2 pc) of ~ 20 X-ray sources around (17 /l 20 m 03 s ,-35°58'20") that can be clearly seen 
in the X-ray image shown in Figure [2] (bottom panel) . This group at the center of the XD cluster 
has the highest surface density in the NGC 6334 complex. 



3.3. Heavily Obscured Cluster XE (HII region NGC 6334 A; FIR NGC 6334 IV) 

XE is the most prominent cluster seen in hard X-rays in the southwest quadrant of the 
NGC 6334 complex. It is only very weakly detected in soft x-rays. XE is spatially associated with 
the radio HII region NGC 6334 A and the far-infrared stellar cluster NGC 6334 IV. This positional 
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agreement confirms that the hard X-ray b and can readily det ect and reso lve heavily embedded 
clust ers as seen, for example, in Cepheus B (IGetman et al. 120061 ) . W3 Main (jFeigelson &; Townslev 
20081 ). and the Rosette Molecular cloud (jWang et alJbood ). 



XE is a rich and concentrated cluster of obscured X-ray sources associated with far-infrared 
source FIR IV and the bipolar radio source A. X-ray E me d ranges from around 2 to an extraordinary 
Emed ~ 5 — 6 keV, corresponding to logiV# ~ 22.3 to 23.8 cm" 2 or Ay ~ 10 to 300 mag. These 
unusually heavily obscured sources in XE include #607, 669 (= FIR-IV 43), 692 (= FIR-IV 60), 
712 and 715. They likely lie within the dense, ~ 1 pc-diameter molec ular clump emitting brightly 
in high-excitation CS rotational emission lines (jKraemer et al.l ll999). As with the group around 



IRN V-l(Sj37T]), we emphasize that young stars are rarely seen with E mec i > 5 keV in other star 
forming regions (see , for example, the Chandra Orion Ultradeep Project E mec i distribution in 
Feigelson etallhoosh . 



Several of these heavily absorbed sources lie near the millimeter sources identified by ISandell 



(|1999h . Sources #737 and 739 lie near MM la, #712 lies near MM lb and MM 2, and #712 lies near 
MM 2. Source #703 with E med = 3.0 keV (corresponding to \ogN H ~ 22.5 cm " 2 and Ay ~ 20) 



coincides with the infrared star IRS 20 associated with infrared nebula IRN IV 2 (Hashimoto et al 



2008). The region around MM-3 is studied by iPersi et al.l (J2009J); only one of their embedded 
infrared-excess stars appears as an X-ray source (#757 = IRS 9). There is no concentration o f 
sources around the ~ 20" shell-like radio HII region NGC 6334 A discussed by lCarral et al.1 ([2002). 



The hard band X-ray cluster XE has a very high stellar surface density and and inferred total 
stellar population of > 500—1000 stars when undetected stars are considered ( §2.4p . The GLIMPSE 
image of the region shows a complex of irregular dense dark clouds superposed on bright emission 
nebulosity. The elongated obscu ring material resemble Infrared Dark Clo uds (IRDCs) noted in 
other star forming regions (e.g. ICarey et al.l Il998l : iRathborne et al.l 120061 ). Some of the heavy 
absorption seen in X-ray stars may be due to these dense obscuring irregular clouds. 



3.4. Extragalactic Blazar RCM B 



The very bright X-ray source #787 coincides with the compact radio continuum source NGC 
6334 B (= RCM B = IGR J17204-3554) which has been widely discussed as a radio-loud active galac- 



tic nucleus lying behing t he NGC 6334 star forming region ([Rodriguez et al.l Il988l : iTrotter et al 



1998 



Bassani et al. 



2005 



Bvkov et al.l 120061 ) . This source is an order of magnitude too bright for 



even the most extreme flaring pre-main sequence stars, and has no stellar counterpart corresponding 
to a cluster OB star. Its ACIS median energy E me d = 4.7 keV (Table 1) indicates heavy absorption 
corresponding to logNfj = 23.3 cm" 2 (Ay ~ 100 mag) consistent with the source lying behind the 
NGC 6334 cloud and other Galactic Plane material. 
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3.5. Lightly Obscured Cluster XF 

NGC 6334 XF is a small and sparse cluster ~0.5 pc across seen best in soft X-rays. It lies 
2 pc (projected) north of the star forming ridge that defines the NGC 6334 complex. Cluster 
membership appears to be dominated by source #821 (K s = 5.71), #846 [K s = 7.36), and #881 
(= CD -35°11482 = IRAS 17170-3539, K s = 7.76, B0.5e). Source #881, with >500 ACIS counts, is 
the brightest X-ray source in the NGC 6634 complex after the 07 star #972 = FIR-III 13 in cluster 
XH (£ |3.7p . The unusually bolometrically bright stellar counterpart of # 821 may be a supergiant 
like HD 319701 and HD 319702 ( pTLl) . If correct, this suggests that the XF cluster is older than 
most of the other clusters with its most massive member evolved off of the main sequence. 

X-ray absorptions of XF members are typically in the range E me d ~ 1.5 — 3 keV corresponding 
to log Njj ~ 21.9 - 22.5 cm -2 and Ay ~ 5 - 20 mag. The Spitzer GLIMPSE image of this region 
(Figure [7]) shows very dark absorption bands characteristic of IRDCs. The GLIMPSE image shows 
bright nebular emission throughout the vicinity (previously labeled NGC 6334 I(NW), see Figured] 
bottom left panel), perhaps with an IRDC superposed. 

3.6. Lightly Obscured Cluster XG 

This cluster in similar to XF, lying off of the molecular ridge and appearing only in the soft 
X-ray image. The cluster is somewhat richer than XF with ~ 25 X-ray sources, implying a total 
population around 250 — 500 members. It is the only x-ray cluster that is located to the south 
of the star forming ridge of NGC 6334. Source #911 is an unusually bright X-ray emitter in the 
cluster. Its stellar counterpart is a heavily reddened star with R — K = 6.6 mag from DENIS 
catalog photometry. Cluster members have a range of absorptions centered around E me d ~ 2 keV, 
corresponding to log Nh ~ 22.0 cm -2 and Ay ~ 5. 

3.7. Heavily Obscured Cluster XH (HII region NGC 6334 D; FIR NGC 6334 II) 

XH has the highest stellar surface density in the hard X-ray image. It is moderately heavily 
absorbed with typical E me( i ~ 2 — 3 keV corresponding to log Nh ~ 22.0 — 22.5 cm~ 2 or Ay ~ 5 — 20 
mag. A few of its sources appear in the soft X-ray band on the southwest side; this might be a 
separate sparse cluster. The X-ray concentration accurately coincides with the far- infrared source 
NGC 6334 II and, displaced ~ 1' east of the X-ray peak, with the large radio HII region NGC 6334 D 
(Figured]). This HII region appears as a bubble- like structure in the GLIMPSE map (Figure [7]) 
and other tracers of heated interstellar material. 



The brightest star in XH seen in both the near-infrared and X-ray bands is source #972 
coincident with the K s = 7.39 star known variously as TICID 73 74-463 from the Tycho Input 
Catalog, FIR-III 13 from the near-infrared photometric survey of IStraw Hvlandl (|l989l ). and 
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IRN III- 1 from the near- infrared polarimetric survey of ([Hashimoto et al.ll2008l ). With 1653 ACIS 
counts, it is the brightest X-ray source in NGC 6334, appearing several arcminutes off-axis of both 
ACIS fields. It was marginally detected as the Einstein Observatory source 2E 1717.1-3548. As it 
lies 2' southwest of the dense core of XH and is only moderately absorbed with E me d = 1.8 keV, 
this star may be better interpreted as a isolated O star rather than as a member of XH. It may be 
the principal ionizing source of the northern portion of the NGC 6334 C HII region. 

Other bright stellar counterparts to X-ray sources in XH include #1291 (= FIR-II-23) and 1312 
with K s ~ 9, and #953 (= FIR-III 10), 1051, 1107, and 1275 (= FIR III 24) with K s ~ 10 - 11. 



3.8. Heavily Obscured Cluster XI (HII regions NGC 6334 E and F; 

FIR NGC 6334 I) 

NGC 6334 XI is a small, sparse cluster seen in the hard X-ray band along the molecular ridge of 
NGC 6334. It corresponds closely to the far-infrared source I (= IRAS 17175-3544) and is centered 
between the two radio continuum HII regions NGC 6334 E and NGC 6334 F (Figured]). The HII 
region NGC 6334 E, centered an arcminute north of the X-ray cluster, has a nearly circular radio 
continuum shell with ~ 30" diameter and a compa ct radio source at the middle associated with a 
very heavily absorbed near-infrared star TPR 161 (jCarral et alj |2002). This massive star does not 
appear in the Chandra image and there is no apparent concentration of lower mass stars within the 
radio HII region. 

Although the XI cluster appears centrally concentrated in the smoothed maps of Figures \M5\ 
examination of the individual source positions in Figure [2] suggests the cluster may be bifurcated 
with subclusters centered at (17 /l 20 m 52 s , -35°45'50°) and (17 ft 20 m 53 s , -35°46'20"). The southern 
X-ray subcluster in cluster XI lies directly on the unusually strong far-infrared sour ce NGC 6334 I 
with a 60 /im flux density ~ 11,000 Jy and a bolometric luminosity ~ 2 x 10 5 L Q ( Sandelll |2000| ) . 
This is the far-infrared source NGC 6334 I and IRAS source 17175-3544. A cluster of near-infrared 



stars extends over 2' (jTapia et al.lll996l ). 



At millimeter wavelengths, the core of NGC 6334 I is resolved into several sources, some 



assoc iated with a cometary ultracompact HII region in radio continuum and masers ([Hunter et al 



2006). This is clearly a concentrated subregion of current massive star formation. X-ray sources 
# 1428, 1429, 1430, 1431, and 1435 lie within a ~ 3" region associated with NGC 6334 I - SMA3 , 
IRS-1E, IRS-I 1 and the bright head of th e rad io continuum structu re. See iHunter et al.l (|200a ). 
Beuther et all (|2007l ) . beuther et al.l (|200J), and lSeifahrt et al.1 (|2008l ) for detailed multiwavelength 
maps of this crowded subregion. Two X-ray sources in this crowded region have counterparts within 
1" at long wavelengths. Source #1428 with E mec [ = 2.8 keV (corresponding to log Nh ~ 22.5 cm" 2 
or Ay ~ 20 mag) coincides with the star IRS-I- 1 (Keck II OSCIR designation, iDe Buizer et al 



2002) or IRS1-W (Very Large Telescope NACO designation, ISeifahrt et al.ll2008l ) with estimated 
spectral type B2. Source #1429 with E me d ~ 3.0 keV coincides with IRS1-SW seen with the 
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NACO camera (jSeifahrt et al.ll2008l ). Several other X-ray sources are present close to this dense 
core, and the locally high density of photons assigned to detector background strongly suggests 
that additional fainter sources would be resolved in a longer X-ray exposure. 



The cometary ultracompact radio HII region NGC 6334F (= G351.42+0.64; ICarral et alJl2002l ) 
lies ~ 10" west of the peak concentration of NGC 6334 I and is associated with the far-infrared 
source IRAS 17175-3544. The moderately strong and absorbed X-ray source #1406 lies within the 
~ 6" radio nebula and possibly associated with the near-infrared star TPR 43. It could be an X-ray 
detection of the exciting star of NGC 6334 F. About 10 unassigned X-ray photons lie within a few 
arcseconds of #1406, representing either the shocked stellar wind or emission from unresolved faint 
low-mass stars. 

On the northern edge of the X-ray cluster XI lies the important protostellar region called 
NGC 6334 I(N). This is a strong extended (sub)millimeter continuum emitter associated with 
M ~ 2000 Mq of de nse molecular gas, molecular outflow and maser emission, but no radio or 
infrared counterparts (jSandelll 120001 . an d references therein ). It has been resolved int o a compact 
cluster of millimeter massive protostars (jHunter et al.ll2006l ). See lBeuther et al.l ((2008) for detailed 
multiwavelength maps of this subregion. Three faint X-ray sources are detected in the ~ 30" 
NGC 6334 I(N) subregion: #1447, 1454 and 1463. Sources #1454 and 1463 exhibit extraordinarily 
high X-ray absorption with E me d ~ 5.5 keV, corresponding to logNjj ~ 23.5 cm -2 and Ay ~ 200 
mag. These are very likely members of the embedded cluster. Source #1454 lies within l"of the 
millimeter source SMA 7 and a faint Spitzer IRAC star, while #1463 may be associated with the 
near-infrared star TPR 216 and/or methanol masers. Source #1447 appears coincident with SMA 6 
but, due to its lower X-ray absorption (E me( i ~ 1 keV), may instead be a foreground star. 



3.9. Lightly Absorbed Cluster XJ 

NGC 6334 XJ is a sparse and small cluster (~ lpc across) seen only in soft x-rays. Its dominant 
star is the counterpart of #1528, 2UCAC 16657515 with V = 14.2 and K s = 9.3. 



3.10. Absent Clusters 

The long-established far-infrared source NGC 6334 III at the center of NGC 6334 is thought 
to be more evolved than sources at the southwest and northeast extremeties. It has no masers 
or outflows but is associated with the large, diffuse bipolar HII region NGC 6334 C. The infrared 
source lies between the lightly-obscured X-ray cluster XG and the heavily obscured cluster XH, but 
there is no concentration of X-ray sources clearly associated with it or the NGC 6334 C HII region. 
We suggest that no stellar cluster is present here. 
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3.11. Unclustered massive stars 



A number of X-ray sources with 8 < K < 11 counterparts lie outside the main cluster con- 
centrations. Several bolometrically luminous stars appear as X-ray sources west of the XA-XD 
clusters. The faint X-ray s ource #2 co incides with HD 319701 (K=6.76), an isolated and lightly 
obscured Bllb supergiant (|Neckelill978l ) lying ~ 10' southwest of the XA-XD complex. The faint 
source #63 similarly lies southwest of the main clusters and appears to be associated with a heavily 
absorbed massive star. Several other sources bright both in the X-ray and K s bands are distributed 
in this region are likely previously unnoticed massive stars (e.g., # 77, 79, 108 and 130). 

The bolometrically brightest of these stars are counterparts of #564 with a K s ~ 8.1 coun- 
terpart lying ~ 2' northeast of XB, #949 with a K s ~ 8.6 counterparts lying ~ 2' north of XG, 
and #972 near cluster XH discussed in ^3.71 Other apparently massive stars distributed across the 
field include # 628, 719, 726, 728, 748, 757, 817, 890, 953 and 1051. These may be massive or 
intermediate-mass stars that have drifted or been ejected from one of the rich clusters. 

Source #1371 coincides with the X-ray and bolometrically bright star HD 319702 {K s ~ 7.39) 
on the eastern side of the NGC 6334 complex, se yeral arcminu tes south of cluster XI. It has been 
spect roscopically classified as a Bllb supergiant (INeckell Il978l ) and a 08 III((f )) giant (IWalborn 
19821 ). Although HD 319701 an d HD 319702 we re identified as the exciting stars of the NGC 6334 
HII region in early work (e.g., ISharplessI Il959l ). it now appears likely that they are older stars 
that have drifted from the main concentrations of heated molecular gas. Other optically bright 
O stars commonly associated with NGC 6334 (HD 156738, 319699 and 319703) lie outside the 
Chandra fields of view. HP 319703 is a bright X-ray source detected with the Einstein Observatory 
(|chlebowski et al.lll989h . 



4. DISCUSSION AND CONCLUSIONS 

This Chandra-based study of the stellar population in the NGC 6334 star forming region must 
be viewed as preliminary in several respects. First, the construction of the sensitive catalog of 
X-ray sources has deficiencies. Some of the faintest sources may be spurious, source properties 
such as X-ray luminosities are not estimated, and inhomogeneities in the survey coverage are not 
quantitatively treated. Second, the exposures are short so that only 5 — 10% of the full Initial 
Mass Function is detected. Third, the interpretation of X-ray photons which are not assigned 
to resolved sources is not addressed here: do they arise from fainter X-ray st ars which must b e 
present in the clusters ( §2.4p . or from diffuse plasma from shocked O star winds (jEzoe et al.ll2006l )? 



Both the incompleteness of the X-ray catalog to faint cluster members, and the interpretation of 
these unassigned photons, require substantially longer Chandra exposures than available for this 
study. Fourth, the association of X-ray sources with near-infrared sources is inadequate. The 
2MASS catalog has insufficient resolution to resolve crowded cluster regions and to obtain reliable 
photometry in the presence of spatially complex nebular emission. Scientific inferences from infrared 
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color-magnitude and color-color diagrams, such as Initial Mass Functions and protoplanetary disk 
fractions, can not be reliably obtained. 

Despite these deficiencies, our survey represents the deepest penetration into the Initial Mass 
Function of the young population over the full NGC 6334 region, and is effective in identifying both 
dozens of luminous OB stars (many of which were not previously identified) and hundreds of lower 
mass members of concentrated stellar clusters. These advantages accrue in large part to the low 
contamination of the X-ray selected population by extraneous sources. The principal results are: 

1. A catalog of 1607 faint X-ray sources is presented, most with subarcsecond positions. Al- 
though only half have 2MASS counterparts, we infer from past studies that ~ 95% of these 
sources are likely cluster members. The remaining sources are extragalactic and stars unre- 
lated to NGC 6334. The X-ray survey suffers far less contamination from Galactic field stars 
or diffuse nebular emission than infrared surveys, and should be nearly unbiased with respect 
to protoplanetary disk emission. In addition, we obtain estimates of line-of-sight absorption 
to each source individually from the median energy of its detected photons. 

2. The X-ray exposures are roughly complete to L x ~ 1 x 10 30 erg s _1 . Based on X-ray studies 
of the stellar populations in the nearby Taurus and Orion clouds, this sensitivity captures the 
top ~ 5 — 10% of the pre-main sequence X-ray luminosity function and most of the cluster 
members with masses above 1 M . The total inferred population of NGC 6334 is around 
20 - 30, 000 stars. 

3. The spatial distribution of X-ray sources is very inhomogeneous and differs greatly when 
lightly obscured sources are compared to heavily obscured sources. From smoothed maps of 
the source distribution, ten X-ray selected clusters are identified, labeled XA to XJ. Seven of 
the clusters are mostly seen in the soft X-ray band and, with out exception, were previously 
unnoticed in visible, infrared or radio studies. The three hard-band X-ray clusters, plus a 
small group behind clusters XC-XD, are closely associated with well-studied luminous far- 
infrared clusters. HII regions are loosely associated with these clusters, but are not generally 
concentrated in the center of the clusters. 

4. The X-ray population includes some of the previously identified OB stars, and detected a 
considerable number of additional stars with K s band magnitudes indicative of new OB 
members. Some of the massive stars are dispersed across the X-ray fields, several parsecs 
from the concentrations of current star formation. At least two appear to be (super)giants. 
This suggests that these massive stars have either drifted or been ejected from their natal 
sites, and that star formation has been ongoing in NGC 6334 for several million years or 
longer. 

5. The X-ray population also includes some extraordinarily obscured sources with X-ray absorp- 
tion up to log Nh — 24.0 cm -2 or Ay ~ 500 mag. It is likely that such high absorptions 
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cannot be produced by int erstellar cloud mater ial, but more likely by localized infalling en- 
velopes around protostars (jGetman et al.l 120061 ) . Several of these sources emit mostly in the 
submillimeter band, an indicator of Class classification. These include well-studied systems 
in NGC 6334 I, NGC 6334 I(N), NGC 6334 IV, and NGC 6334 V. If some of these X-ray 
sources are truly massive Class protostars known from submillimeter studies; they may be 
the youngest stars ever detected in the X-ray band. The only comparably embedded systems 
reported in the X-ray band are the lower-mass Class 0-1 sources IRS 7 in Coronal Australis 
cloud with logN H ^ 23.5 c m" 2 and IRAS 21391+5802 with logiVff ~ 24.0 cm" 2 in the 
triggered cloudlet IC 1396N (|Hamaguchi et al.1 120051 : iGetman et al.l 120071 ) . We note, how- 
ever, that it is difficult with the available data to distinguish emission from the submillimeter 
source itself from closely associated lower mass pre-main sequence stars. For example, the 
Becklin-Neugebauer Object in the OMC-1 core behind the Orion Nebula is a very faint X-ray 
source lying only ~ 1" (~ 400 AU) from a much more X-ray luminous lower mass companion 
(jGrosso et alJbood ). 



It is clear that there is great promise in further study of the stellar population of NGC 6334 
based the results obtained here. First, the existing 40 ks Chandra exposures are photon-starved, 
with most detections having <10 counts and the vast majority of lower mass stars undetected. 
A > 10- fold increase in exposure would alleviate this problem and provide an empirical basis for 
interpreting the unresolved X-ray emission as stellar or diffuse. Second, wide-field, near- and mid- 
infrared images with sub-arcsecond resolution are urgently needed for comparison with the X-ray 
source catalog. The challenge here is less reaching high sensitivity, as the counterparts to the X-ray 
sources will mostly have K s < 16, as resolving stars from the bright, spatially complex nebular 
emission. Third, near-infrared spectroscopy of X-ray selected stars would provide spectral types and 
direct estimates of stellar masses. The highest priority for spectra are the 8 < K < 11 counterparts 
which are likely new OB members of the cluster. 
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Fig. 1. — Mosaic of two ACIS-I images of the Chandra of NGC 6334 with 1607 source extraction 
regions shown in yellow. 




Fig. 2. — Two expanded views of crowded portions of the Chandra fields. Top: Far-infrared regons 
NGC 6334 IV and V and the new clusters XA-XD. Bottom: Far-infrared NGC 6334 I and II regions 
and the new cluster XJ. See Figure [5] for the cluster identifications. 
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Fig. 3. — Near-infrared properties of candidate 2MASS counterparts of the Chandra X-ray sources. 
Left: Color-color diagram for counterparts having high-quality photometry. Red triangles show 
stars with colors consistent with reddened photospheres (Class III stars), and green circles show 
stars with ET s -band excesses associated with protoplanetary disks (Class II stars). Dotted lines 
show reddening vecto rs marked every Ay = 5 mag. The locus of ma in sequence stars is from 
Bessell fc Brett (119881 ') , classical T Tauri stars from Meyer et al.1 (119971 ) , and HAeBe stars from 
Lada &i Adamsl (119921 ). Right: J vs. J-H col or-magnitude diag ram of these stars. 1 Myr and 2 Myr 
isochrones for pre- main sequence stars from lsiess et al. ( 2000 ) are shown in addition to the ZAMS. 
Note that only a small fraction of the 1607 Chandra sources are shown in these diagrams due to 
limitations of the 2MASS catalog. 
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Fig. 4. — The spatial distribution of the 1607 X-ray sources shown as a stellar surface density map 
after smoothing with a a = 45" Gaussian kernel. 
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Slellor Surface Density (Log ^stars/qrcmin') 




Fig. 5. — The smoothed spatial distribution of soft X-ray sources (first panel) and hard X-ray 
sources (second panel). The sample is divided at E me d = 2.5 keV which approximately corresponds 
to column density log Njj ~ 22.3 cm -2 and visual absorption Ay ~ 10 mag. The ten X-ray clusters 
are outlined and labeled. 
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Fig. 6. — The X-r ay source ma p supe rposed on long-wavelengths maps. Top: an 18 cm radio 
contin uum map by ISarma et al.1 (120001 ). Bottom: a 71 //m far- infrared map by lLoughran et al. 
(jl98a ). The left and right panels show the soft and hard X-ray stellar surface density maps, 
respectively. 
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Fig. 7. — Spitzer GLIMPSE map with the X-ray fields and cluster superposed. Here the 3.6/um, 
5.6/im, and 8//m intensities are represented as blue, green and red, respectively. [7] 



Table 1. Chandra Catalog of X-ray Sources in NGC 6334 
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Note. — Tablc^is available in its entirety in the electronic edition of the Journal. Only a portion of the table around NGC 6334 I(N) is shown here. 
Col. (1): X-ray catalog sequence number, sorted by RA. 
Col. (2): IAU designation. 

Cols. (3) and (4): Right ascension and declination for epoch (J2000.0). 

Col. (5): Estimated standard deviation of the random component of the position error, \j a % + a y ■ The single-axis position errors, a x and cr y , arc estimated from the single-axis 
standard deviations of the PSF inside the extraction region and the number of counts extracted. 
Col. (6): Off-axis angle. 

Cols. (7) and (8): Net counts extracted in the total energy band (0.5—8 keV); average of the upper and lower lcr errors on col. (7). 
Col. (9): Background counts expected in the source extraction region (total band). 
Col. (10): Net counts extracted in the hard energy band (2—8 keV). 

Col. (11): Fraction of the PSF (at 1.5 keV) enclosed within the extraction region. A reduced PSF fraction (significantly below 90%) may indicate that the source is in a crowded 
region. 

Col. (12): Log of the probability that the source arises from fluctuations in the local background based on Poisson statistics. 
Col. (13): Background-corrected median photon energy (total band). 

Col. (14): K s magnitude of 2MASS Point Source Catalog stars within 1 of the X-ray position. ":" indicates photometric confusion by a nearby star or nebular emission, or 
photometric accuracy worse than ±0.2 mag. ">" indicates source was detected in the J or H bands but not the K s band. 

Col. (15): Possible counterparts within Z" of the X-ray source in the SIMBAD database of published objects. FIR-V sources are from [St raw fc Hvlandl i 19891) , TPR. stars arc from 
Tapia ct al ] i 19961 1. and other designations arc discussed in the text. 
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Table 2. X-ray Star Clusters in NGC 6334 
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Note. — 

Col. (1): X-ray Cluster designation 

Col. (2) and (3): Right ascension and declination for epoch (J2000.0) 

Col. (4)and (5): Indicator of cluster richness in soft (col. 4) and hard (col. 5) band: Y = Yes, 
N = No, W = Weak 

Col. (6): The spatial extent of the x-ray cluster, rounded to the nearest arcminute. 

Col. (7): X-ray source count in the X-ray cluster. 

Cols. (8): IRAS point source catalo g counterpart, nr = n ear. 

Cols. (9): Far-infrared counterpart I Loughran et al.l 1986 ) 

Cols. (10): Radio HII region counterpart IjRodriguez et al.Ml982h ■ 



